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We evaluate the effect of screening by bound electron in 7 Be(p,7) 8 B reaction, where 7 Be target 
contains bound electron, in the framework of the adiabatic representation of the three particle 
problem. A comparison with two other approximations (united atom and folding) is presented. A 
good agreement between the "united atom" approximation and the exact solution is found. We 
also discuss the screening corrections induced by two K-shell electrons on a 7 Be target. The bound 
electron screening effect consequences for 7 Be and 8 B solar neutrino fluxes are discussed. 



I. INTRODUCTION 



In recent years, an increasing attention has been devoted to an accurate estimation of electron screening effect for 
. nuclear fusion reactions in stellar plasma and for the interactions of low-energy ion beams with atomic or molecular 
targets in laboratory experiments (see refs. and references therein). 

In this Letter we present the first quantum mechanical calculation of screening effect by bound electron in 



7 Bc + p — > 8 B + 7 (1) 



nuclear fusion from the pp-cycle in the Sun. Contribution of this reaction into the the total luminosity of the Sun is 
negligible small, but it is directly related to the long-standing "Solar Neutrino Problem" , - one of the most intriguing 
issue in the present-day neutrino astrophysics. Standard physics cannot explain an 37 Ar production rate in the 
£C) • Chlorine experiment smaller than that expected from the solar 8 B neutrino flux measured by both Kamiokande and 
(with better statistic) Super Kamiokande. GALLEX and SAGE experiments also indicate beryllium neutrino deficit 
(see the discussion in ref. ||). One of the most elegant solutions to the solar neutrino anomaly is resonant neutrino 
flavor conversion in the sun, that is the so-called MSW effect j p~0| . It requires an extension of the minimal standard 
electroweak theory: neutrino masses and neutrino mixing. These neutrino oscillation parameters are determined in a 
CIh' wav that can bridge between the predictions of the standard solar models and the solar neutrino observations. Thus, 
even in the framework of the standard solar model within a hypothesis of neutrino oscillation (and MSW effect), it 
' is apparently needed more precise calculations of nuclear fusion rates in the sun, because they can significantly affect 
\ the neutrino oscillation parameters determination. 

A careful study of electron screening effect on nuclear fusion rates becomes particularly actual in view of expected 
high accuracy neutrino flux measurements by a number of new large detectors (Super Kamiokande, SNO). The 
interpretation of forthcoming data requires relevant precise calculations of solar neutrino fluxes and neutrino energy 
y\ ' spectrum. 

Usually, the effects of surrounding plasma on the nuclear fusion are treated in electrostatic screening approximation. 
This approximation, being classical or quantum, correctly reflects the major properties of a process only for high 
relative velocity of the colliding nuclei, when electron density in the vicinity of the fusing nuclei remains almost 
unchanged during the collision. In the case, when relative velocity of the nuclei is much smaller or comparable 
with the electron one (and this is the case at solar conditions), the electron density changes following any relative 
configuration of the nuclei, and the electrons have an impact on a kinetic energy shift of the nuclei. It is therefore 
natural to consider the phenomenon within the framework of the adiabatic approximation, which comes from the well 
known Born-Oppenheimer (BO) approximation. 

The BO approximation allows one to treat nuclear motion independently from the electron coordinate, within a 
new effective potential which depends on the internuclear distance. Since, the considered nuclear velocities are smaller 
than the boun electron one, corrections to the BO approximation are expected to be negligible. Obviously, the fusing 
nuclei are from the continuum energy spectrum. An accurate treatment within the adiabatic approximation of the 
screening effect by electrons from continuum spectrum requires an additional research but the case of bound electrons 
presents no special problem (see, for example ref. |pd|). 

As it was argued by A. Dar, G. Shaviv, and N. Shaviv ^,^2|, the commonly accepted Debye-Hiickel theory is not 
quite adequate for evaluating the screening effect in not- very-dense stars, like the Sun. There is, also, an experimental 
evidence that this theory does not provide correct answer for the screening || . Actually, this fact is of no importance 
when the screening due to the plasma electrons is by itself rather small. But it is not the case for the low-lying bound 
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electrons which do screen the electric charge of nuclei much effectively, and moreover, the screening effect drastically 
increases when energy of the fusing ions decreases. The electron screening can have dramatic effects in very dense 
stellar cores. 

At low and moderate energies, the fusion cross section of "bare" charged nuclei colliding with the relative momenta 
p in the center-of-mass frame is expressed as (see ref. Q): 

fTp\ -27TJ7 /r>\ 

^ ' = ~~E~ ' ^ ' 

where S(E) is the so-called astrophysical factor which incorporates all nuclear features of the process, E is the collision 
energy of the nuclei, 77 = M ZiZ 2 / '(m e a p) is the usual Coulomb parameter, m e and M are the electron and reduced 
nuclear masses, respectively, and ao is the hydrogen Bohr radius. The exponential factor originates from the Coulomb 
wave function of the internuclear motion ipg(R) at R = 0. 

The screened cross section a s (E) differs by the enhancement factor 

where i^e{R) is the wave function of the internuclear motion which accounts for the bound electron. 

We evaluate the effect of electron screening of 7 Be nucleus by one bound electron in reaction in the framework 
of the adiabatic approximation for three particle problem. This calculation is compared with two relevant approx- 
imations, "united atom" (UA) and folding approximations which, as we will demonstrate below, give respectively 
upper and lower estimates for the screening effect. In the framework of the UA approximation we estimate also the 
screening effect for 7 Be nucleus with two K-shell electrons. 



II. METHOD OF CALCULATION 



We treat the Coulomb problem for three particles in the framework of the adiabatic representation, that consists 
as follows. The three particle Schrodinger equation solution ^(r, R) is represented as 

00 

*(r,R) = ^^(R)0 fc (r;R). 

Here, R is the internuclear radius-vector, and r is the electron radius-vector from the center-of-mass of nuclei. The 
two-center eigenf unctions ^>fc(r;R) are derived from the Schrodinger equation for three particles, when two nuclei are 
fixed on a distance R, and the eigenfunctions have a dependence on R as a parameter. The electron energy eigenvalue 
U n i m (R) also depends on the internuclear distance. The nuclei with electric charges Z\ and Z 2 interact with the 
effective potential 

V nlm (R) = + U nlm (R), (4) 

K 

and ^fe(R) are the eigenfunctions of nuclear motion in the potential We shall write the parameters in atomic 
units, except that we shall display the units. Thus, the energy unit is m e e 4 /ft 2 r* 27.21 eV, the lenght unit is ao, and 
the electric charge unit is e. 

Our approximation consists in using only one two-center eigenfunction corresponding to the ground state of the 
system. Obviously, the main screening effect is reached at zero orbital moment of the colliding nuclei. Thus, the 
fusing nuclei are considered in S-wave. There are some arguments for this approximation. At first, the high energy 
states corrections (at fixed R) are of the order of \J m e /M. Then, excited energy levels correspond to the less energy 
of the united nucleus. It leads to their exponentially small contribution into the nuclear fusion rate in comparison 
with the ground state of the system. At least, only the ground state energy of the electron in the field of two nuclei 
( 7 Be and p in our case), labeled as ISa therm (see, for example ref. |H) has the correct asymptotic behaviour (see 
below). Thus, the total three particle wave function reads: 

*(r,R) =^(RMr;R). 

The tabulated values of U n im(R) |@ are used for our purposes. The ground state eigenvalue energy U(R) (= Uqoq{R)), 
with zero quantum numbers, has the correct asymptotic behaviour: the electron eigenvalue energy reaches the energy 
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of the united ion as R approaches zero: U(0) = (Z\ + Z 2 ) 2 /2, and U(oo) = Zf/2 (that is the energy of the isolated 
ion eZi). In fig. the values for the ground state are plotted. 

We solved numerically the Schrodinger equation for two scattering nuclei with Z\ = 4 and Z2 = 1 at center-of-mass 
kinetic energy E within the potential (|J) with zero quantum numbers n = I = m = 0. The numerical solution was 
obtained on the mesh on R for R = and R = i? ma x in Numerov's scheme. We varied integration step h and R max in 
order to ensure that final result is not changed substantially. Also for checking purposes we reproduced the value of 
^(0)| 2 substituting U(R) — 0. The values of the squared wave function of relative motion of the two nuclei |-0£;(O)| 2 
at 0.1 keV < E < 100 keV are calculated and compared with |V^(0)| 2 = 2m]/(e 2m > - 1). 

The UA approximation can serve as an upper limit to the screening corrections induced by bound electron. The 
essence of this approximation is replacement of the kinetic energy E by the shifted energy E + AE, where AE is 
the energy difference of the bound electron between the final and initial states. For one bound electron this energy 
difference is given by AE = {Z x + Z 2 ) 2 /2 - Z 2 /2. 

Thus, the UA approximation changes VB(0)| 2 to \iPe+ae(®)\ 2 ■ Since the classical turning point for the ionic barrier 
penetration is at much smaller radius than the Bohr radius (a = a /Zi) at solar energies, the nuclear fusion takes 
place, in the main, within the united ion state. The classical turning point is given by: 

ZiZ 2 /0.4keV\ 
Rb = ^* a [— E—)> 

therefore, at solar energies, the UA approximation applicability is well satisfied. Moreover, at kinetic energies above 
few keV, U(R) is not changed substantially at < R < Rb, and the UA approximation is expected to be very close to 
the exact solution. The UA approximation replaces decreasing U(R) (see in fig. [l]) by the constant AE, hence, this 
approximation serves as an upper estimate. 

The folding approximation can be used as a lower limit to the screening corrections. This approximation embodies 
the static wave function of the bound electron, which does not depend on the internuclear distance. The total three 
particle wave function in this approximation is presented as a product: 

* / (r,R)=^ / (R)0 / (r) > 



where the wave function for the electron bounded on the nucleus with electric charge Z\ reads: (f>f(r) — \J Z\j-Ke^ Zir , 
and ipf(R) is the wave function of the two nuclei. Averaging the electron coordinates, ipf(R) can be derived from the 
Schrodinger equation with the potential: 

Vf(R) = ^-Z 2 -( i^-^ - Zl e- 2 «A . (5) 



R \ R 

The zero point of the colliding nuclei kinetic energy is the electron eigenvalue energy in the field of Z\ nucleus. 
Obviously, a rigid electron wave function, which has no dependence on R, can be used only in the case of instantaneous 
impact of the nuclei. Thus, the area of the folding approximation applicability is the faster velocity of the nuclei relative 
to the electron one. Attractive part of the potential Vf reaches the value: Z\Z 2 as R — ► 0. Since Z\Z 2 < (Zi + Z 2 ) 2 /2 
that is the value of the electron energy in united ion, the folding approximation gives a lower estimate for the 
enhancement factor. 



III. SUMMARY AND DISCUSSION 

The enhancement factor (||) is plotted in fig. ^ for all three approximations. As it was expected, the UA approxi- 
mation always overestimate the exact solution, while the folding approximation underestimates it. Nevertheless, it is 
easy to see that simple UA prescription gives very close values to the exact solution at kinetic energies above 2 keV. 
Therefore, the latter can be used not only as a qualitative, but as a good quantitative approximation to the electron 
screening by bound electrons. 

The electron screening is dramatic at very low kinetic energies of the nuclei. However, in a plasma, most of the 
nuclear fusions come at the Gamow peak energy, that is defined by the strong dependence of the nuclear cross section 
on energy (||) and the fast decrease of the exponential particle distribution. This energy is given by: 

E = 1.22(Z 1 Z 2 T 6 ) 2 / 3 (Af/A/ p ) 1 / 3 keV, 

where Tq = T/10 6 K, M p is the proton mass. In the solar interior at r e s/RQ — 0.06, where the 7 Be and 8 B neutrino 
production reaches its maximum p6[ , the plasma parameters are T 6 rj 14.7, the electron density n e rs 7.7/qq, and 
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the Gamow peak energy in reaction (Q) is about 18 keV. Then 7(-Eo) = 1-1, that is, there is 10% of an enhancement 
by one bound electron in the bohron production rate. Simple computations within the UA approximation give the 
screening effect as: 



^AE/kT 



(6) 



Then, 10% of an enhancement by one bound electron could be easily reproduced just inserting numbers into the 
formula (||). One can apply this formula also for a 7 Be nucleus with two bound electrons. Then, 

A£ = (x? + X?) - (X? c + X2°) = 227.98 eV. 

Here, xf = 340.2 eV and %f = 259.4 eV are, respectively, the fifth and forth ionization potential of the 8 B atom 
and xf e = 217.72 eV and %2 6 = 159.9 eV are, respectively, the forth and third ionization potential of the 7 Be atom. 
Thus, two bound electrons enhancement factor is given by j{Eq) = 1.196, i.e. roughly 20%. 

Using the Saha equation Iben, Kalata and Schwartz jlTj calculated the probabilities fx and fi that one or two 
K-shell electrons are associated with any given Be nucleus. The calculations were perfomed under the assumption 
of pure Coulomb electron-ion forces, neglecting all excited states and screening. The probabilites found are 



/i = A 



I - A i n. 25 A 2 exp ( ^ 



where 



2-KmJtT 



3/2 



exp 



Xi_ 
kT 



Here k is the Boltzmann's constant, and A x = \i — X2 = 63.8 eV. Inserting numbers one can obtain: fx = 30%, 
and f 2 = 3%. 

Using the calculated abundances of 7 Be ions, one can estimate the thermal averaged screening effect induced by 
both one and two bound electrons on a 7 Be nucleus: 



(7) - 1 = 0.30 x 0.1 + 0.03 x 0.2 « 0.04. 

In the Standard Solar Model (SSM) the electron capture rate by 7 Be nucleus is taken to be about 1000 times faster 
than the proton capture rate |L6| ]. Therefore, a small change in 8 B production rate does not affect significantly the 
7 Be neutrino flux, although it makes a proportional change in 8 B neutrino flux. Thus, the electron screening by bound 
electrons has the prompt consequences on bohron neutrino flux. 

The electron screening by plasma electrons from the continuum spectrum is expected to contribute significantly to 
the total enhancement factor, since it is proportional to n e , and thus it has to be taken into account in the exact 
prediction of neutrino flux change. 

In summary, in the solar interior K-shell bound electrons enhance Be(p,7) 8 B rate and increase 8 B neutrino pro- 
duction rate by of about 4%. Therefore, bound electron screening has an effect on the solar 8 B neutrinos, and acts 
with the opposite effect to the berrylium neutrinos. The main essence of the electron screening in nuclear fusions is 
the change in electron density on a nucleus during the collision of the nuclei. This effect can be treated only in a 
dynamical calculation like the present three particle calculation or the discussed UA approximation. 

We thank J. N. Bahcall, A. V. Gruzinov and V. A. Naumov for usefull discussions. 



[1] A. V. Gruzinov and J. N. Bahcall, |astro-ph/9801028 . 

[2] T. D. Shoppa, M. Jeng, S. E. Koonin, K. Langanke and R. Seki, Nucl. Phys. A605 (1996) 387; T. D. Shoppa, S. E. Koonin, 

K. Langanke and R. Seki, Phys. Rev. C 48 (1993) 837. 
[3] E. E. Salpeter and H. M. Van Horn, Ap. J. 155 (1969) 183. 
[4] H. E. Mitler, Ap. J. 212 (1977) 513. 



4 



[5] C. Carraro, A. Schaefer and S. E. Koonin, Ap. J. 331 (1988) 331. 

[6] L. S. Brown and R. F. Sawyer, Rev. Mod. Phys. 69 (1997) 411. 

[7] C. W. Johnson, E. Kolbe, S. E. Koonin and K. Langanke, Ap. J. 392 (1992) 320. 

[8] N. Shaviv and G. Shaviv, Ap. J. 468 (1996) 433. 

[9] S. DegPInnocenti et al., in Proc. of the Workshop on Trends in Astroparticle Physics, Stockholm, Sweden, 1994, ed. 
L. Bergstrom et al. [Nucl. Phys. B (Proc. Suppl.) 43 (1995) 66]. 
[10] L. Wolfenstein, Phys. Rev. D 17 (1978) 2369; D 20 (1979) 2634; S. P. Mikheyev and A. Yu. Smirnov, Yad. Fiz. 42 (1985) 

1441 [Sov. J. Nucl. Phys. 42 (1985) 913]; Nuovo Cimento 9 C (1986) 17. 
[11] L. Bracci, G. Innocenti, V. S. Melezhik, G. Mezzorani and P. Pasini, Phys. Lett. A 153 (1991) 456. 
[12] A. Dar and G. Shaviv, Ap. J. 468 (1996) 933. 
[13] Astrophysical formulae (Springer, Berlin, 1974). 

[14] S. S. Gershtein and V. D. Krivchenkov, Sov. Phys. JETP (1961) 1491. 

[15] L. I. Ponomarev and T. P. Puzynina, Dubna preprint R4-3175 (1967), in Russian. 

[16] J. N. Bahcall and M. H. Pinsonneault, Rev. Mod. Phys. 67 (1995) 781. 

[17] I. Iben, K. Kalata and J. Schwartz, Ap. J. 150 (1967) 1001. 



5 




2 4 6 8 10 12 14 16 18 

R 

FIG. 1. Electron energy in the field of the two nuclei 7 Be and p 




FIG. 2. Enhancement of nuclear fusion rate due to the electron screening 
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